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ABSTRACT 

Outflows from young stellar objects have been identified as a possible source of turbulence in molec¬ 
ular clouds. To investigate the relationship between outflows, cloud dynamics and turbulence, we 
compare the kinematics of the molecular gas associated with NGC 1333, traced in 13 CO(1-0), with the 
distribution of young stellar objects (YSOs) within. We find a velocity dispersion of ~ 1 — 1.6km/s in 
13 CO that does not significantly vary across the cloud, and is uncorrelated with the number of nearby 
young stellar outflows identified from optical and submillimeter observations. However, from velocity 
channel maps we identify about 20 depressions in the 13 CO intensity of scales > 0.1 — 0.2pc and 
velocity widths 1 — 3km/s. The depressions exhibit limb brightened rims in both individual velocity 
channel maps and position velocity diagrams, suggesting that they are slowly expanding cavities. We 
interpret these depressions to be remnants of past YSO outflow activity: If these cavities are presently 
empty, they would fill in on time scales of ~ 10 6 yr. This can exceed the lifetime of a YSO outflow 
phase, or the transit time of the central star through the cavity, explaining the absence of any clear 
correlation between the cavities and YSO outflows. We find that the momentum and energy deposi¬ 
tion associated with the expansion of the cavities is sufficient to power the turbulence in the cloud. 

In this way we conclude that the cavities are an important intermediate step between the conversion 
of YSO outflow energy and momentum into cloud turbulent motions. 

Subject headings: ISM: kinematics and dynamics — ISM: individual objects: NGC 1333 — ISM: 
clouds — ISM: bubbles — ISM: jets and outflows — ISM: molecules 


1. INTRODUCTION 

Energetic outflows exert a strong effect on their 
parent molecular clou ds (for a recent review see 
iBallv fc Reipurthll20f)T ). Measurements of the total ki¬ 
netic energy present in young stellar outflows imply that 
outflows and winds associated with young stellar objects 
contain sufficient kinetic energy to excite a significant 
fraction of the supersonic turbulence present i n molecular 
clouds and/or unbind and disper s e the cloud lllallv et_alj 
H196Ij_|BaIl y_fcReiDurthl 1200 ll iKnee fc Sanded! 120001: 
Warm et al .11 1996th 

The molecular cloud associated with the reflection 
nebula NGC 1333 provides a good setting to study 
the relation between young stellar outflows and their 
influence on parent clouds because it is an active 
site of low and intermediate mass star formation. 
NGC 1333 contains two young star cluste r s identified 
in ne ar-infrared studies llAsnin et ahl ll 9941 1 1 .ad a et ahl 
dU), many Herbig-Haro ( HH) objects and associ¬ 
ated molecular outflows (e. g., IBallv et "allIl996el lAsninl 
120031 IKnee fc Sandellll2000D as well as younger, embed¬ 
ded IRAS sources llSande ll fc Knedl2000l Pfennings et all 
110871 iRodrignez. et al .11 19901) . 

Previous studies have suggested that the outflows 
present in NGC 1333 could have perturbed the associated 
molecular cloud. Based on his study of the star clusters, 
lAsninl J20f)l suggested that the large number of active 
molecular outflows in the southern region of NGC 1333 

Electronic address: aquillen@pas.rochester.edu 


could provide a mechanism for exciting cloud turbulence 
and causing subsequent st ar formation. Based on a study 
of the cloud in 12 COf3-2b (Knee fc Sandell l200(tl found 
that molecular outflows observed in the cloud could pro¬ 
vide enough kinetic energy to accelerate the entire cloud 
by a few km/s. Dust ridges and shells seen in submil¬ 
limeter continuum can be associated with outflows imply¬ 
ing that the cloud itself has been considera bly m odified 
and, perhaps, disrupted bv the outflows llSand ell fc Knee 


200(1 iWarin et all 1 19961: iLefloch et~allll998t IBallv et al 


1996a). 


Many authors have argued that outflows could stim¬ 
ulate star formation or be responsible for turbulence 
in molecular clouds. These arguments usually rely 
on estimates of the energy present in outflows com¬ 


pared to that associated with turbulence in the clouds 
llBallv & Reipurthll200ll l. Studies of individual objects, 
however, express a more complicated picture. In partic¬ 
ular the explicit nature of the coupling between cloud 
material and outflows could be better explored. Molec¬ 
ular outflows are identified through millimeter radio 
emission from molecules such as CO extending beyond 
the more collim ated emission seen in HH objec ts in 
the visible (e.g., iBence et~ mum is; et, al . 1999/) ■ In 
giant outflows, CO emission associated with the out¬ 
flow ca n extend to a few areminutes from the central 
sourc e llChernin fc Massnlll 9961 iCernicharo Hr. Reinurthl 
[T99fl . Unified models suggest that it is the narrow 
bipolar jet with, perhaps, a wide-angle wind compo¬ 
nent, that entrains molecular material which is then 
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seen as the CO molecular outflow dRag ajg_Ca^rijll993: 
IChernin et all 11994 IDelamarter et all I2000t ILee et al l 

l2002t lArce fc Sargent! 12004 iGardiner et aLlboOrih . We 

note that high angular resolution observational studies 
have revealed evidence for wide angle winds in a ddition 
to collimated jets (e.g. . lArce fc Goodmanll2002aj) . 

In 12 CO lines such as the J =1-0 transition the op¬ 
tical depth increases rapidly near the velocity centroid 
of a molecular cloud. However, the optical depth de¬ 
creases with increasing velocity offset and so 12 CO can be 
used to trace lower density material asso ciated with out¬ 
flows (e.g., the survey of NGC 1333 bv lKnee fc SandeU 
mj. However to trace the structure at higher densi¬ 
ties, a mo re op tically thi n tra cer such as 13 CO or C 18 0 
is needed. IR.idge etaT. ( 120031) have carried out a survey 
of nearby molecular clouds in these two emission lines, 
providing a database with which to investigate the mech¬ 
anism and timescale for the dissipation of gas around 
individual stars and clusters. NGC 1333 h as been pre¬ 
viousl y surveyed in the same molecules by IWarin et a,l.l 
ifTiiol . The overall structure of the warm gas in the 
cloud is expected to be well traced by the C 18 0 emis¬ 
sion, with depletion or freeze-out of CO only occurring 
in small dense pockets (iBergin et alii2002 ). 

Outflows and winds associated with young stellar ob¬ 
jects produce high velocity gas, but also can evac¬ 
uate regions in their host molecular cloud. Cav¬ 
ities associated with outflows have been seen pre¬ 
viously, pa rticularly in low opacity tracers such as 
13 CO (e.g.. IWarin et alJll996HArce fe Goodmanll2002bl: 

iLefloch et all 1 1 99 8t IWelch et all 120001) . In the Circi- 
nus outflow. iBallv et alJ (119991) described dusty filaments 
confining lpc sized cavities that were also seen in their 
13 CO and C 18 0 maps. They suggested th ese cavities 
were fossil remnants of previous outflows. Snell et all 
iTmsnl) described a 15 krn/s velocity, pc sized do uble 
lobed cavity associated with L1551. iReinurth et alJ 
iTiosl) described an evacuated pc long chimney associ¬ 
ated with HH-310. With higher resolution CO inter¬ 
ferometry, ILee et afl (120021) described complex outflow 
morphology such as shell-like and multi polar st ru ctures, 
multip le cavities, and asymmetric lobes. lArce fc Sargentl 
: 200 1) showed that the outflow of L1228 had evacuated a 
conical shaped cavity, and had morphology and kinemat¬ 
ics consistent with entrainment of gas by bot h coll i mate d 
and wide angle outflow components. IWelch etld . (j200 (f) 
showed that a 0.05pc expanding shell in Taurus could 
have been caused by outflows d riven b y the binar y XZ 
Taurus. In the NGC 1333 cloud. IWarin et all ( 1996 ) sug¬ 
gest that the SVS13, IRAS 2 and IRAS 4 outflows have 
created the cavity that l i es bet ween SVS13 and HH12. 
However iKnee fc Sandell 1 2000 ) found that the molecu¬ 
lar outflow associated with IRAS 4 did not overlap with 
this cavity. 

In a few cases, theoretical models for expanding spher¬ 
ical cavities i.e. wind blown bubb les, have bee n applied 
to cavities opened by outflows. iKoo fc McKee) 1 992 a: 
applied their wind blown bubble model to the HH 7- 
11 region, located in the NGC 1333 c loud. In the con¬ 
text of a slow wind, IKoo fc McKeel (Il992al) suggested 
that cavities opened in molecular clouds by outflows were 
consistent with their estimated mechanical luminosities. 
iRaga et all ( 120041) modeled the limb -b rightened cavity 
associated with HH 46/47 (iNoriega-Crespo et ali 1200 4) 


as a bow shock driven by a perfectly collimated jet. 

Most observational studies have suggested that exist¬ 
ing HH objects and outflows could be responsible for cre- 


ating cavities (e.g., 
1992al Warm et al 

i\nee 6 

199<f 

l sanaeii zuuu: ivoo <v; ivicivee 

Lcfloch et al. 199^ Lg^et al. 

2002 

Arce & Goodman 

20021/f Norieg^i-Crespo et al. 


2004, IRaga et all [20041) 7 and have identified the most 


prominent outflows evident in the winds of the CO chan¬ 
nel maps or seen as HH objects with the most prominent 
cavities. Since the duration of the strong outflow period 
of a protostar may be short compared with the life of 
a young cluster, the history of the coupling between the 
outflow and molecular cloud must be considered. Indeed, 
some studies have also suggested that cavities in molec¬ 
ular clouds could be tracers of previously active outflows 
(e.g. . IBallv et a,U fl999). In the vicinity of S VS13 in NGC 
1333 , both possibilities have been suggested (IWarin et alJ 
H996tlLefloch et, aLlll998l) . IWarin et al.l (119961) suggested 
that the star clus ter north of S VS13 was associated with 
a cavity, whereas fLefloch et alJ ( 1998 ) identified two cav¬ 
ities south of SVS13 and associated one them with a 
cone shaped outflow emanating from the infrared source 
SVS13. 


To investigate the coupling between outflows and the 
dynamics of molecular clouds, we compare the spatial 
distribution of outflows in NGC 1333 with the properties 
of its molecular cloud. Specifically we compare the lower 
opt ical depth 1 3 CO millimeter velocity cubes obtained 
hv lRidge et al (120031) with near infrared sources and the 
locations of outflows identified by previous surveys. In 
addition we compare structures seen in the molecular gas 
distribution to the locations of young stellar objects. 

We adopt a distance to NGC 1333 of 212 pc based 
on the distance estimat ed to t he sta r BD+30°549 that 
illuminates the nebula dCerni ima. This distance is 
nearer than the estimate 318 pc that is based on on the 
Hipparcos p arallactic measuremen ts of the Perseus OB2 
association (Ide Zeeuw et al . nu. At the distance of 
212 pc, 1’ corresponds to 0.062 pc. 


2. MEAN VELOCITY FIELD AND VELOCITY DISPERSION 
OF THE MOLECULAR GAS 

If outflows associated with young stellar objects are 
responsible for turbulence excited in the cloud, then we 
might expect regions that contain large numbers of out¬ 
flows to reside in a more turbulent molecular medium. 
In this section we compare the integrated pro perties of 
the molecular gas, as traced in 13 CO by (iRidge et alJ 
12003 ) with the distribution of previously identified out¬ 
flow sources. 

Of the NG C 1333 v elocity cubes presented and dis¬ 
cussed bv IR.idge et all ( 2003 ) the highest signal to noise 
data set for NGC 1333 is the 13 CO(1-0) data cube ob¬ 
tained at the Five College Radio Astronomy Observa¬ 
tory (FCRAO). Because of its higher signal to noise we 
primarily use this velocity set for our comparison. The 
angular resolution of these data is 47" and the spectral 
resolution corresponds to a velocity spacing between each 
channel of 0.133 km/s. For NGC 1333, 47" corresponds 
to ~ .05pc. 

In Figure H we show the mean velocity measured from 
the 13 CO(1-0) data cube as contours, overlaid on the 
2MASS Ks band image shown as grayscale. In Figure 
[3 we show the velocity dispersion measured from the 
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13 CO(1-0) data cube, overlaid on the integrated inten¬ 
sity of the velocity cube. The mean velocity, < v >, is 
the intensity weighted mean velocity at each position of 
the map measured from the entire velocity cube and is 
computed as < v >= JA V ili/Y2i h where and V{ are 
the velocity and intensities of the individual pixels. The 
velocity dispersion, cr, is the square root of the intensity 
weighted square of the velocity subtracted by the mean; 

* 2 = £(«<-<« >) 2 V *■ . 

We might expect an increase in the velocity dispersion 
of the molecular cloud in region that contain many out¬ 
flow sources or many embedded young objects. However, 
from inspection of FigureQwe find that the velocity dis¬ 
persion as measured from 13 CO does not significantly 
vary across the molecular cloud. The velocity dispersion 
ranges between 1.0 and 1.6 km/s across the cloud. This 
velocity range is well within the blue-shifted wings seen 
in 12 CO that are 10 km/s below the cloud velocity mean 
from ou tflows associated with the most prominent HH 

Hydrostatic equilibrium models for clouds predict re- 
latio ns between velocity d ispersion, density and radius 
(e.g., iMvers fc Fulled ir992'l. However, Fig. 0 shows no 
strong relation between integrated intensity and veloc¬ 
ity dispersion, as would be expected from hydrostatic 
cloud models. There is little significant structure in the 
constant dispersion contours, though the dispersion is 
somewhat higher near the lower 13 CO intensity region to 
the north of SVS-13. The somewhat higher velocity dis¬ 
persion in this depression or cavity might be explained 
with a scenario in which a cavity is evacuated by winds 
driven from sources within the cavity. Two star clus- 
ters have be en identified in the near-infrared imaging by 
[Aspin et"aH 1 1994 1. One cluster is located in the vicinity 
of this cavity and another is located near the depression 
in the integrated 13 CO about 7’ to the north of SVS13. 
These clusters can be seen in the K-band stellar distri¬ 
bution in Fig. □ The outflow region associated with 
SVS-13 (the HH7-11 region) does not have a larger ve¬ 
locity dispersion than other regions of the cloud. Like¬ 
wise the HH4, HH2 and HH12 regions a s sociat ed with 
outflows seen in 12 CO bv lKnee fc Sandell 1200(111 also do 
not stand out in the velocity dispersion maps. Significant 
structure correlated with positions of known outflows in 
the cloud random motions as traced by the velocity dis¬ 
persion are not observed in the the lower opacity tracer 
13 CO. 

Outflows previously identified in 12 CO as high velocity 
wings (with blue and red-shifted emiss ion greater than 
a few km/s from the cloud mean) by l|Knee fc Sandelll 
MM) are not evident in the high velocity wings in the 
13 CO or C 18 0 data of llRidge et alJ I 200 .lll . At veloci¬ 
ties below 5 km/s and above 10 km/s little emission is 
detecte d in the individual 13 C O channel maps (see Fig¬ 
ure 0. iKnee fc Sandelll 1200011 found outflows associated 
with HH 6 , HH12, HH7-11, IRAS 2 and IRAS 4; their CO 
maps did not extend all the way to IRAS 1. Surveys in 
12 CO detect outflows w hereas those in 13 CO often do not 
(e.g.. lWarin et al.lllll996 l. This is most likely because the 
12 CO emission is more sensitive to optically thin and low 
density regions, but less cap able of penetrating de nser 
higher opacity regions (e.g., lArce fc Sargent! 12004 s ). A 
deeper, more sensitive velocity cube in 13 CO might re¬ 
veal additional structure near the individual outflows, 


particularly at velocities more than 3 km/s above or be¬ 
low the mean cloud velocity. Outflows traced in the high 
velocity wings of molecules such as 12 CO are associated 
with a low density medium that is faint and so difficult 
to detect in 13 CO. 

Because we see little high velocity gas in 13 CO and 
no large spatial variations in the velocity dispersion, we 
find little evidence for turbulence in the denser regions 
of the molecular cloud that has been directly induced by 
the high concentration of high velocity outflows discov¬ 
ered at other wavelengths. We conclude that the bulk of 
the cloud, which we expect is traced in 13 CO, does not 
contain much high velocity gas associated with active 
outflows. 

In contrast to the velocity dispersion map, there is 
quite a bit structure in the mean velocity field shown 
in Figure ^ The mean velocity field exhibits a gradient 
with velocity increasing from the south to the north of 
the cloud by about a km/s. There is also a jump in ve¬ 
locity of ~ 0.3 km/s separating the southern region in 
NGC 1333S from the rest of the cluster. 

3. CAVITIES 

In the previous section we found that the integrated 
properties of the 13 CO velocity cube were not strongly 
influenced by outflows. In this section we inspect individ¬ 
ual channel maps and position velocity diagrams from the 
13 CO velocity cube to better investigate the possibility 
that outflows have affected the structure of the molecular 
cloud. 

Instead of high velocity gas associated with outflows, 
we see shells and cavities or depressions in the gas distri¬ 
bution in individual channel maps; see Figure^ Cavities 
in the 13 CO integrated intensity have been described by 
previ ous studies (e.g., IWarin et al.l |199(I iLefloch et all 
11998 1. Because these depressions in intensity are seen 
in more than one neighboring velocity channel, they are 
likely to be real and not a result of noise in the data. 
Different cavities appear at different velocities, suggest¬ 
ing that they are often not directly related to one an¬ 
other. Some cavities can disappear at velocities higher or 
lower than the central value where the depression is most 
prominent, although CO emission can remain strong. 
This suggests that the cavities can lie within the cloud, 
and are not exclusively a result of filamentary structure 
near the cloud edges. Previous studies comparing maps 
of line emission in C 13 0, 13 CO, NH 3 and CS have found 
that depressions in the line intensity in 13 CO and C 18 0 
are primarily due to lower gas density rather than de- 
pletion of CO, temperature variation s or high opacity 
ijWarin et a,l ill 996t ILefloch et alJI1998ll . From inspection 
of the channel maps, we find that the most prominent 
cavities are also evident in the lower opacity and lower 
signal to noise C 18 0 velocity cube of NGC 1333 (for ad¬ 
ditional discussion on morphology in C 18 0 compared to 
that seen in 13 CO see section 3.1 below). 

The velocity gradient across the NGC 1333 molecular 
cloud is fairly smooth or shallow (Fig. EJ, so evacuated 
or lower density regions stand out as depressions in indi¬ 
vidual channel maps. The rims of many of these evacu¬ 
ated regions are oval in the low angular resolution 13 CO 
images (see Fig. QJ. If molecular gas has been evacu¬ 
ated from a central region then we would expect that 
the resulting cavity would appear limb brightened. This 
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appears to be true for a number of cavities evident in the 
individual channel maps (e.g., Cl, C4). 

In Figure 0] we have identified and labeled a number of 
cavities. Their estimated central positions and velocities 
are listed in Table 1 in order of increasing central velocity. 
These cavities correspond to depressions in channel maps 
that we have also identified as depressions in position 
velocity diagrams. To list a cavity, we required that they 
be visible (by eye) as depressions in the intensity seen 
in both channel maps and position velocity diagrams, 
however we did not require a specific contrast level in 
the intensity. In some cases the intensity of the edge of 
the cavity is only ~ 50% that of the center in a particular 
channel, in other cases the intensity inside the cavity is a 
few times lower than that of its rim (e.g., CIO). We also 
required that the depression be surrounded on at least 
three sides by higher intensity emission in a channel map. 

In Fig. 0 we show position velocity diagrams extracted 
from the 13 CO velocity cube in narrow horizontal strips. 
The width of the strip corresponds to one pixel in the 
velocity cube or 25" which is less than half the beam 
width. The x —axes of these plots is the RA(J2000) and 
the y— axis, the velocity in km/s. Figure 0 shows posi¬ 
tion velocity diagrams extracted for a number of different 
horizontal (East-West strips) each separated by 2' in dec¬ 
lination. The DEC(J2000) of each strip is shown in the 
upper right hand side of each subplot. Figure El shows 
a similar set of position velocity diagrams, however each 
plot is extracted along a narrow strip oriented vertically. 
Consequently the x-axes are DEC(J2000). The separa¬ 
tion in RA between strips is 2'. The RA(J2000) of each 
strip is shown on the upper right hand side of each sub¬ 
plot. 

3.1. Cavity Properties 

The cavities that we have identified in Figures El El 
and El have typical widths of 2-4' corresponding to 0.1- 
0.2pc. The density contrast suggested by the difference 
between the emission intensity at the center of the cavity 
compared to that in the edge brightened rims is a factor 
of 2-4. Cavities are evident in channel maps and are also 
seen as depressions in the position velocity diagrams. For 
many cavities, at the cavity center emission can be seen 
red-shifted from the mean central velocity and also blue- 
shifted from this mean. If the cavity is expanding, the 
red-shifted emission we would interpret to be behind the 
cavity moving away from us and the blue-shifted would 
correspond to material in front of the cavity moving to¬ 
ward us. The width in velocity space of these cavities is 
approximately 1 km/s except for extreme cavities such 
as CIO which is seen in the channel maps between 6 and 
9 km/s. Some cavities are bounded by both red-shifted 
and blue-shifted emission (e.g., Cll) others appear to 
have one open end (e.g., CIO). 

Many cavities are connected to or associated with other 
cavities. In other words they are near in position and 
velocity to other cavities. For example, Cll, C12, and 
C13 together may form one long tubular region oriented 
north-south in the cloud. C9 and C14 are a smoothly 
connected triangular region, and C7 and C4 could be 
connected. C15 and C5 could be a connected tubular 
region running south-east. In the southern region of 
the cloud, the interconnected cavities seem to run north- 
south, whereas on the eastern side of the cloud, cavities 


tend to be oriented east-west. The orientation of the 
cavities would be consistent with a scenario where the 
cavities were associated with previous or relic outflows 
driven from sources near the center of the cloud. The 
center of the cloud contains two clusters of young stars 
which are evident in near-infrared images (see Figure 0. 
These sources could have driven active outflows a few 
hundred thousand to about a million years ago (lAspinl 
12001 . 

Cavities that have red-shifted and blue-shifted rims 
(e.g, C13, Cll, C9, C3, C9, C8, Cl, C7, C16) are most 
easily interpreted in terms of expansion. The radial ex¬ 
pansion rate is slow, half to 1 km/s typically, however 
this is high enough that is above the sound speed for 
the cold temperatures typical of molecular clouds ~ 0.3 
km/s, though it is similar to the cloud’s velocity disper¬ 
sion. A few of the cavities (e.g., C5, CIO and C19) have 
larger radial expansion velocities (1.2-2.5 km/s), above 
the velocity of turbulent motions in the cloud. The red- 
shifted and blue-shifted rims would be consistent with 
a cavity that is slowly expanding in the cloud. Because 
emission can be brighter in the red-shifted or blue-shifted 
rims, compared to nearby material (e.g., C9, C16), the 
cavities are likely to be expanding rather than boundaries 
of regions of the cloud which have only been evacuated 
by outflows. 

An expanding spherical shell would appear as an oval 
symmetrical about the spatial axis in position velocity 
diagrams. If the position velocity diagram is extracted 
along the major axis of an expanding ovoid shell, the cav¬ 
ity would appear as a skewed oval (one oriented along a 
diagonal line) in the position velocity diagram, with one 
side r ed-sh ifted_and the other the other side blue-shifted 
(e.g., LWeIdLefoal,| |2I)a. If the position velocity dia¬ 
gram is extracted along a different axis, the cavity would 
also appear as an oval in the position velocity diagram 
but with central position and velocity dependent upon 
the location of the extracted position velocity diagram. 
Since we were searching the entire cloud for cavities, we 
have extracted in Figures El El position velocity diagrams 
in north-south and east-west strips in the cloud. If non- 
spherical (elongated or ovate) cavities are oriented at dif¬ 
ferent positions with respect to the line of sight and the 
angle that the position velocity diagrams are extracted, 
we would we expect to see skewed ovals in the position 
velocity diagrams as well as ovals with positions that vary 
with positi on in neigh boring position velocity diagrams 
(e.g., see [Welch et al,l l2000h . 

In the position velocity diagrams we see examples of 
cavities which could be interpreted as expanding non- 
spherical or ovate shells at different orientations. For ex¬ 
ample, CIO appears to be oriented nearly along the line 
of sight. It is a small cavity that spans 3 km/s and has al¬ 
most no skew in the position velocity diagrams. In other 
words, the cavity central position is not strongly depen¬ 
dent upon the velocity channel, and the central velocity 
does not strongly depend on position on the sky. How¬ 
ever, C18 and C19 together may be an example of a cav¬ 
ity that is skewed in the position velocity diagram, with 
the western side red-shifted and the eastern side blue- 
shifted. The connected cavities (e.g., Cll, C12, C13) 
could be examples of cavities which change location in 
both velocity and position in different position velocity 
diagrams. It is likely that we are seeing expanding ovoid 
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cavities oriented at different angles with respect to the 
viewer as well as those oriented along the line of sight, 
such as CIO. CIO is seen at almost all velocities (except 
near 5 krn/s) suggesting that it passes through the entire 
cloud. Other cavities, such as C12 and C15 are seen only 
at higher or lower velocities, suggesting that they pass 
through only one side of the cloud. 

A few cavities exhibit high velocity wings (of order 
2 km/s from the mean which is above the velocity dis¬ 
persion of the cloud). CIO is extraordinarily broad, ex¬ 
tending through a number of channel maps. The cav¬ 
ity has a width of about 3km/s. C5 could be con¬ 
nected to blue-shifted high velocity gas (at 4 km/s), 
see the panel at RA=3h29ml9s(2000) in Fig. El The 
evacuated triangular shape in the position velocity dia¬ 
gram suggests that a conical shaped cavity could have 
been opened. This shape has b een seen in previous 
works. For example, lArce "fc Sargent ( 1201) 1 1 suggested 
that IRAS20582+7724 exhibited a cavity which had been 
opened by a wide angle wind. C19 is also a triangular 
shaped cavity in the position velocity diagram (see the 
panel at DEC= +31°18 , 04"(2000) in Fig. EJ and is as¬ 
sociated with red-shifted gas 2km/s above the mean of 
the cloud. Because they are associated with higher ve¬ 
locity emission on their rims, CIO, C5 and C19, could be 
examples of younger cavities. 

The C 18 O(l-0) velocity cube has lower signal to noise 
than the 13 CO(1-0) cube, however the C 18 O(l-0) line 
also has lower optical depth. Features detected in 13 CO 
in the outer more tenuous regions are more difficult to 
detect in the C 18 0 velocity cube. However, cavities iden¬ 
tified in the denser regions of the cloud should be visible 
in both velocity cubes. In Figure0we show position ve¬ 
locity plots of the C 18 0 emission extracted in the same 
regions as those of Figure El which displays the 13 CO(l- 
0) emission. From a comparison between Figure Eland 
0we see that cavities such as C6, C7, Cll, C12, C13, 
C14 and C16 which have higher intensity rims are seen 
in both velocity cubes. The cavities which are not de¬ 
tected in the C 18 0 velocity cube lie in the more diffuse 
regions of the cloud. Particularly prominent in the C 18 0 
position velocity plot is the cavity C16. The good cor¬ 
respondence between the morphology seen in these two 
figures confirms that the cavities identified in 13 CO in 
the denser regions of the cloud are not due to noise in 
the data, and are not artifacts introduced by variations 
in optical depth. 

3.2. The relation between cavities and infrared sources 

In this section we compare the location of the cavi¬ 
ties seen in individual channel maps to the positions of 
sources identified from previous surveys. A comparison 
between the location of cavities and the near-infrared 
Ks-band 2MASS images showing T-Tauri stars, young 
stellar clusters and background sources is shown in Fig¬ 
ure |H1 We have plotted our cavity locations along with 
the positions of sources identified from infrared imag¬ 
ing surveys, llAsDinl [200^JLaRaetalJ lT99fit IStrom et al.l 
fl97l . HH objects IIReiPurtbl 119991) and molecular out¬ 
flows studied bv lKnee fc Sandelll (1200011 . The positions 
of the B stars BD+30547 and BD+30549, and infrared 
source SVS13 are also shown on this figure. 

While it was tempting to associated every cavity with a 
source, no matter how faint, in practice some cavities did 


not contain sources likely to be responsible for them. For 
example, cavities on the southern side of the cloud such 
as Cll, C12 have no HH or 2MASS counterparts. HH 
objects and embedded young stell ar objects (as seen from 
the submillimeter observations of ISandell fc Kneell200(1lf 
appear to reside within the cloud and are not located in 
these cavities. The outflows bright in shocked CO and 
molecular hydrogen and those identified as HH objects, 
tend to be located along regions of bright 13 CO emission 
rather located in cavities. For example, IRAS 1 is located 
at a peak in the 13 CO emission in most channel maps 
rather than a depression in any channel map. 

The lack of correspondence between the cavities and 
molecular outflows, HH objects, submillimeter and IRAS 
sources, suggests that the cavities are not directly asso¬ 
ciated with the youngest stellar objects. We now con¬ 
sider the possibility that the cavities are related to older 
class II sources which should be visible in the 2MASS Ks 
band images. Stars in the stellar clusters evident in the 
2MAS S im ages or the infrared imaging surveys of lAsoinl 
l(200l and lLada et all (^996) could be associated with 
ca vities near th e cen ter of the cloud. T he stars studied 
by lAspinl (I2003D and lLada et ahl (11996f) (identified from 
near-infrared imaging) have ages approximately a 1 mil¬ 
lion years old. Because of the large number of the stars 
seen in Ks band on the northern side of the cloud, it 
would be possible to match each cavity in this region 
with a pre-main sequence object. However the lack of 
near-infrared objects on the southern and eastern regions 
of the cloud makes it impossible to associate each cav¬ 
ity with a nearby pre-main sequence object. While the 
lack of correspondence between younger sources (embed¬ 
ded with outflows) and cavities implies that the cavities 
could be relics of previous activity, the lack of correspon¬ 
dence between cavities and older sources (see in the near- 
infrared images) suggests that cavities could be created 
distant from the current location of the source originally 
responsible for driving it. 

If the cavities are relics, then the stars which could 
have caused them may have drifted from their location 
of birth. T-Tauri stars hav e been estimat ed to be approx¬ 
imately a million years old lAspinl ll2003i) . A star moving 
at 1 km/s (approximately the velocity dispersion of the 
cloud) would be able to move 0.1 pc in 10 5 years. This 
is far enough that it could have moved across a cavity. If 
the cavities are relics of outflows, their source stars could 
have moved away from the location where they drove out¬ 
flows. For example, the embedded class II source IRAS 
16316-1540 lies at the edge of a cavi ty (associate d with 
RN091) that it could have caused lLee et al . lt200l . Cav¬ 
ities could also be relics of long jets (greater than 0.5 
pc long) that have punctured tubular holes through the 
molecular cloud. Thus cavities could be formed distant 
from the source originally responsible for driving the jet. 

In the previous section we pointed out three cavities, 
(CIO, C5 and C19), which have higher expansion veloci¬ 
ties than the others. Since these three cavities are prob¬ 
ably younger, it may be easier to find candidate objects 
responsible for causing them. C5 and CIO are located in 
a region that has been covered by most previous surveys 
whereas C19 is not, making it difficult to search for a 
candidate object responsible for this cavity. A variety of 
types of candidate sources are located in the vicinity of 
C5. CIO contains pre-main sequence stars but no IRAS 
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source, CO outflow, submillimeter source, or HH object 
(see Figure 0. Since CIO is probably pointed toward 
us, we expect the candidate driving source to be in the 
vicinity of the cavity. This suggests that CIO was driven 
by one of the nearby pre-main sequence stars and not a 
younger source with an active outflow. 

3.3. Observational summary 

Before we discuss explanations and scenarios account¬ 
ing for the cavities we have identified in NGC 1333 in 
13 CO we summarize their observational properties. 

1) The molecular cloud of NGC 1333 is full of limb- 
brightened shells and cavities which are evident in both 
channel maps and position velocity slices in the 13 CO 
velocity cube. We have identified 22 cavities in an ap¬ 
proximately 1 x 1 pc region on the sky. 

2) We can estimate the volume filling factor of these 
cavities in the cloud. We have 22 cavities with approx¬ 
imate radii of O.lpc in a volume approximately (lpc) 3 . 
We estimate that approximately 10% by volume of the 
molecular cloud is comprised of cavities. Lower rim ex¬ 
pansion velocity, lower density contrast and smaller cav¬ 
ities would be more difficult to detect in the 13 CO ve¬ 
locity cube. Since we cannot have identified every cavity 
in this cloud, we can regard the estimated filling factor 
as a lower limit. We conclude that cavities permeate the 
molecular cloud. 

3) Limb-brightening is seen in both position space 
(making shells) as well as velocity space, suggesting that 
cavities can be expanding. 

4) Cavities have expansion velocity widths of ~ 1 km/s 
and spatial widths of 0.1-0.2 pc. The intensity contrast 
between that at the rim and that in the center is not high 
~2-4, as estimated from a comparison between the emis¬ 
sion intensity in the cavity center at the cavity’s central 
velocity and the emission temperature at its edge. If the 
edges of the cavities are turbulent then the edge would 
cause emission at a range of velocities. Consequently the 
true density contrast between the density at the rim and 
that inside the cavity could be larger than our measured 
factor. 

5) Cavities appear to be non-spherical. They could be 
ovate or cylindrical. Examples of elongated cavities at 
different orientations are seen in the velocity cube. 

6 ) We see three cavities with higher velocity 13 CO 
emission (above 2 km/s from the cloud mean), suggesting 
that some cavities could be younger than others. There 
could be an evolutionary sequence. 

7) Some cavities appear to be connected to others. 
Tubular cavities are oriented north-south on the south¬ 
ern part of the cloud and east-west on the eastern side. 
This suggests that these cavities could have been driven 
by previously active outflows located near the center of 
the cloud where there are million year old star clusters. 

8 ) Cavities are not directly or obviously associated with 
outflows previously identified as HH objects or seen in 
shocked molecular hydrogen and CO in Spitzer 4.5/xm 
images (Porras, private communication). Instead these 
young stellar objects seem to be associated with dense 
molecular regions or filaments that are emitting (rather 
than deficient) in 13 CO. 

Cavity expansion velocities (~ lkm/s) are observed 
to be above the sound speed of the cloud, though they 
are similar to the size of the cloud’s velocity dispersion. 


For gas traveling at the sound speed of a fraction of a 
km/s typical of a cold molecular cloud, it would take a 
few times 10 5 years for a 0.2pc diameter cavity to fill 
in if the cavity were not expanding and assuming that 
currently there is no momentum source inside the cavity 
adding to the expansion. Turbulent motions in the cloud 
could lower this timescale by a factor of a few. However, 
it might take an additional few times 10 5 years longer to 
fill in the cavity taking into account the need to overcome 
the momentum in the walls implied from the estimated 
expansion velocities. Based on the age of the cluster, 
estimated from the age of the oldest cluster members 
or a few million years, the cavity lifetime is short with 
respect to the age of the entire system. This suggests 
that the cavities have formed fairly recently, even though 
they do not seem to be directly associated with existing 
outflows. 

Our failure to match cavities with existing outflows 
suggests that the cavities identified here are primarily 
relics and associated with past outflow activity, as previ¬ 
ously suggested bv lBallv et all ( 1999 1. Whereas existing 
outflows don’t seem to be strongly perturbing the bulk 
of the cloud, the cavities we identify here permeate the 
cloud. It is possible that we failed to match individual 
outflows to the cavities identified here in part because be¬ 
cause of the low angular resolution of t he 13 CO veloci ty 
cube in comparison to other studies (e.g,. IArce fc Sargentl 
l200l . The 13 CO velocity cube used here is only sensi¬ 
tive to cavities larger than ~ 0. 05 pc since the spatial 
resolution of these data are 47” ijRidge et ai.ll2003h . cor¬ 
responding to ~ 0.05 pc at the distance of NGC 1333. 
We expect that cavities with diameter smaller than O.lpc 
would not only be harder to detect in the velocity cube 
studied here but might also be shorter lived, making 
them rarer. In short, smaller cavities could also be asso¬ 
ciated with molecular outflows and HH objects, but we 
may not be able to detect them in the 13 CO cube studied 
here if they are too small. 

It is possible that the cavities we have identified here 
are are result of random turbulence that is excited by 
energy sources external to the cluster. For example, sim¬ 
ulations of supersonic compressible MHD turbulence in 
a sheering disk can show depressi ons in synthetic chan- 
nel maps (e.g. see Figure 2b of iPichardo et al.l 1200(11 
that mimic the signature of an expanding cavity. Dif¬ 
ferentiating between organized structures such as shells, 
cavities and filaments and those caused by random tur¬ 
bulent motions is a pervasive problem i n the study of the 
ISM. For example lElmegreen fc ScaP (I2004H at the end 
of their recent review of interstellar turbulence ask ’’Why 
is the power spectrum of density structure a power law 
when direct observation shows the ISM to be a collection 
of shells, bubbles, comets, spiral wave shocks and other 
pressurized structures spanning a wide range of scales?” 

To differentiate between structures generated from 
sources within the cloud and those arising from turbulent 
motions (with energy source assumed to be external to 
the cloud), we can attempt to search for differences be¬ 
tween the structures seen in the velocity cube of NGC 
1333 and those exhibited by simulations of turbulent 
media. We qualitatively compare Figure IH to the syn¬ 
thetic channel maps shown in figure 2b bv lPichardo et all 
1200 (1 . Depressions in the synthetic channel maps tend 
to have redshifted or blueshifted rims, but not both (as 
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do Cll, C13 and C16). Depressions in the gas density 
which span a range of velocities larger than the veloc¬ 
ity dispersion of the cloud (such as C5, CIO, and C19) 
appear to be absent in the simulation. This compari¬ 
son suggests that conical shaped evacuated regions might 
be more easily confused with structures created by tur¬ 
bulence, however ovate or spherical or higher expansion 
velocity cavities are likely to have been driven by inter¬ 
nal energy sources. This comparison suggests that some 
(but not necessarily all) of our candidate cavities have 
been driven by energy sources internal to the cloud. But 
if high velocity cavities are likely to have been driven 
by internal sources then we infer that there should also 
be lower velocity structures which correspond to older 
or weaker expanding cavities. Likewise, if cavities with 
both redshifted and blue shifted rims are likely to have 
been driven by internal sources, then we infer that there 
should also be cavities located at the edge of the cloud 
which h ave prominent ri ms only on one side. The simu¬ 
lation of IPichardo et alJ (. 2000 1 illustrates general prop¬ 
erties of supersonic compressible MHD turbulence driven 
by energy at large scales, however the scale of the sim¬ 
ulation (500pc) is larger than that of our cloud (lpc) 
and so the simulation is not an ideal physical analogy to 
NGC 1333. It may in future be possible to develop better 
tools to differentiate between the scenarios for producing 
structures observed in molecular clouds. 


4. ENERGETICS OF CAVITIES 


To investigate the energetics of driving cavities we con¬ 
sider wind blown bubble models (e.g., see the review 
by iFrankl 119991. When the bubble is fully radiative 
we can consider the m omentum-conserving case (e.g., 
:Kq<> fc McKee ll992al'nlK We assume that the outflow 
responsible for the cavity has ceased, consistent with the 
interpretation of these cavities as relics. Therefore, we 
wish to relate the size and expansion velocity of the cav¬ 
ity to the total momentum imparted to the ambient cloud 
from the outflow. By dimensional analysis in the spher¬ 
ical case 

^ - Rh~ l (1) 

P 


where P w is the total momentum output during the life¬ 
time of the outflow; P w = M w V w t w , where M w is the 
mass flux in the wind at a mean velocity V w , during a 
litetime t w . Here R is the radius of the cavity at time 
t and p is the cloud density. This equation is similar to 
the Sedov solution for an explosion and is appropriate if 
the hydrodynamics can be described by scale free func¬ 
tions. Alore_detailedonedimensional spherical models 
(e.g.. lKoo fc McKeelll992albl) support the type of order 
of magnitude scale-free estimate illustrated here. In more 
convenient units we find 


P w ~ 0.2Mokm s 


_i / n{H 2 ) 

\ 10 4 cm -3 


R 


O.lpc 


2 x 10 5 yr 

( 2 ) 

where we have adopted a number density of molecular 
hydrogen n(H 2 ) ~ 10 4 cm -3 typical of the NGC 1333 
cloud. Note this estimate is very sensitive to the radius 
of the cavity. If the cavity is one half the size, then the 
momentum requirements drop by more than an order of 
magnitude. One can also consider the role of the cavity 
age in shaping the cloud. If we assume that all cavities 


are produced by outflows of identical momentum out¬ 
puts, then cavities that differed by a factor of 2 in size 
would differ by a factor of 16 in age. 

For Equation ^ we have assumed a spherical model. 

We can consider a cylindrical cavity with the scaling 

= A 3 * -1 (3) 

pi 

where P Wl ± is the total momentum exerted in the direc¬ 
tion perpendicular to a jet along the jet axis a length l , 
and P w .±/l is the total momentum exerted perpendicu¬ 
lar to the jet axis per unit length. Here l refers to the 
length of the cylinder, which could be half the width of 
the cloud if the driving source is located at the center of 
the cloud and the jet has drilled right through the cloud. 

Using the above scaling relation we would estimate 

P. ± ~ 0.8M s km s -> (-5_ V ( -L-) ( —E-) 

\10 4 cm -3 / \ O.lpc/ \0.4pc/ \2 x 10 5 yr/ 

(4) 

where we have used parameters similar to those in Equa¬ 
tion [3 

Both of the above scaling relations imply that the cav¬ 
ity expands at a velocity v exp = 4? ~ R/t. Different 
scaling factors (1/3 or 1/4) would be inserted in front of 
this relation, depending on the power of R in relations 
Hand0 If the cavities in NGC 1333’s molecular cloud 
are ~ 2 x 10 5 years old we would then predict an expan¬ 
sion velocity close to but somewhat smaller than that 
observed. Consequently for we adopt t ~ 2 x 10 5 years 
since the outflow ceased for the age of the 0.1 pc sized 
cavities in NGC 1333. 

We expect that the ratio of momentum exerted per¬ 
pendicular to the jet axis to that exerted along the jet 
axis would depend on a number of quantities. For a col¬ 
limated jet that with a constant thrust, we expect this 
ratio would primarily depend upon the Mach number of 
the jet. However, this ratio could also depend on the ra¬ 
tio of momentum originating from a disk wind compared 
to that originating from a collimated jet. Young stel¬ 
lar objects probably generate outflows that decay with 
time. Consequently the bow shock would slow down as 
a function of time. The width of the cavity compared to 
its length could increase with time. The cavity’s shape, 
expansion rate at its edges and density contrast (inside 
compared to outside) should be sensitive to the history 
of the outflow responsible for driving it. 

A wind-blown cavity model for these cavities can be 
used to predict observable cavity properties. Using 
Equation [5] we can substitute for the expansion veloc¬ 
ity v exp ~ R/t finding 

Pw ~ pR Vexp ( 5 ) 

where we have placed the observables on the right hand 
’ side and left the total momentum required on the left 
hand side. The observables are the density p, the cavity 
size, R , and the expansion velocity v exp . Using Equation 
Eland making the same substitution, we find 

Pw,A. ^ plR Vexp ■ (6) 

From these relations we see that a cavity’s environment 
could affect its dynamic properties. A similar sized cavity 
in a cloud of higher density should be older than a similar 




















size one in a lower density region and so should have a 
lower expansion velocity. Because it is easier to detect 
large fast cavities in a velocity cube, we suspect that 
cavities will be easiest to identify at the edges of clouds. 
At cloud edges, the lower densities will allow the cavities 
to grow to a large size while maintaining an expansion 
velocity below the channel spacing or that of turbulence 
in the cloud. This does appear to be the case in NGC 
1333, see Figure^ where it is particularly easy to identify 
large cavities near the cloud edge (e.g., C19). 

We also expect scaling relations between cavity observ¬ 
ables. Older cavities should be larger, have lower expan¬ 
sion velocities, and lower density contrasts. We hope that 
future studies which aim to match wind-blown bubble 
models to young cavities (with the highest velocities) can 
better constrain the total momentum requirements as 
well as the nature of the momentum injection; e.g. wide 
angle wind or jet, and momentum injection as a func¬ 
tion o f time llCunningham et aLll2005bt iThorndike et all 

HHI). 

4.1. Comparison of momentum required to open cavities 
to that in observed outflows 

We compare the momentum required to open the 
cavity (estimated from Equation |2l to that estimated 
from young stellar outflows. Young stellar objects (class 
0 and 1 sources) have outflows lasting ~ 10 5 years. 
Momentum fluxes, estimated from observations of en¬ 
trained material, range from 0.1 to 1000 Mmkm/s for low 
and high mass ob jects respectively fHidiei_e^dJ]20(Ki 
iFukui et, alJll993D . In NGC 1333 itself. IKi!ee"faSandeill 
i20odb measured a total momentum, due to about 6 
mostly class 0 sources, of ~ lOM 0 km/s in swept-up high 
velocity gas (here high velocity means > a few km/s from 
the mean). The momentum imparted by the individ¬ 
ual outflows ranged from 0.2 — 2M 0 km/s. We find that 
these total momentum values are similar to those esti¬ 
mated fro m Equ a tionsEl andlllusinK our scaling relations. 
iKnee fc Sandelll ( 20000 measured energies of ~ 10 43 - 
10 44 erg in each outflow which had typical velocities of 
lOkm/s based on red-shifted and blue-shifted 12 CO emis¬ 
sion. 

The momentum traced in the 12 CO bv lKnee fc SandelH 
lil correspond to sources that are younger than the 
estimated ages of the cavities. As the YSOs age we ex¬ 
pect that more momentum will be put into the molecular 
cloud. The rough correspondence between the momen¬ 
tum required to open the cavity and that observed in 
entrained higher velocity 12 CO suggests that a substan¬ 
tial fraction of the momentum from outflows observed in 
12 CO is required to form the cavities we have identified 
in the 13 CO velocity cube. It may be that only a fraction 
of the momentum from outflows is actually ejected from 
the cloud. Thus we conclude that a degree of momen¬ 
tum coupling exists between the outflows and the cloud. 
The difference in number of outflows identified in 12 CO 
(handful) and the number of cavities we have identified 
here (22) implies that cavities can be detected over a 
timescale that is a few times longer that that for out¬ 
flows in 12 CO. 

It is likely that there is only a moderately short win¬ 
dow in time when a cavity would be identified in the 
13 CO channel maps. Cavities with radii smaller than 
0.05pc would be difficult to detect with the 47” beam 


of the 13 CO data. Larger and older cavities, would be 
more slowly expanding and would have lower density con¬ 
trasts. This would make them more difficult to identify 
in the channel maps. Based on Equation [3 we estimate 
that a cavity four times larger could be moving only half 
as fast. Thus we expect a cavity with a 0.4pc diameter 
would have a shell velocity of only 0.2 km/s at which 
time it would become confused with turbulence in the 
cloud and so difficult to identify in an individual channel 
map. (Note the channel maps here have a comparable 
velocity resolution of 0.133 km/s.) Above, we have es¬ 
timated an age of ~ 2 x 10 5 years for a 0.1 pc radius 
cavity. Equation 0 implies that this cavity would have 
about one fourth the expansion velocity when it is twice 
as big, at which time it would be about eight times older 
(using Equation QJ. From this we estimate that an ex¬ 
panding cavity would be detectable in the channel maps 
during a timescale of about a million years. This esti¬ 
mated timescale compared to that for outflows (a hun¬ 
dred thousand years) could account for the difference in 
number of outflows identified in 12 CO and the number 
of cavities we have identified in 13 CO. 

4.2. Excitation of Turbulence 

Previous studies have suggested that outflows from 
young stellar objects are sufficiently energetic or contain 
sufficient momentum that they could account for a sig¬ 
nificant fraction of turbulence in molecular clouds. Here 
we have introduced an intermediate stage, the formation 
of a cavity which is caused by an extinct or previously 
active outflow. In this subsection we compare the total 
momentum required to account for the observed cavities 
to that present in the turbulent cloud. 

The momentum and energy injection rate from 
outflows in NGC 1333 has been estimated to be 
10~ 4 Mf^ km s~ 1 yr~ 1 and O.IL ^ from the outflows seen 
in 12 CO llKnee fa SandelllEiOOOH . In the previous sections 
we found that the momentum required to form the ob¬ 
served cavities is ~ lM^ km/s, similar to that of a typical 
outflow seen in 12 CO bv IKnee' fe SandeH (l2000f) . There 
are a few times as many cavities (22) as outflows (~ 6) 
seen in 12 CO. The number of observed cavities is consis¬ 
tent with the longer lifetime we have estimated for these 
cavities or a few hundred thousand years. Consequently 
the total momentum input into the cloud via the ex¬ 
pansion of cavities is similar to that estimated from the 
observed outflows in 12 CO or ~ 1O _4 M 0 km s _1 yr _1 . 
The total energy flux consist ent with the cavit i es is also 
similar to that estimated by IKnee fc Sandelll 112000H or 
~ O.1L 0 = 4 x 10 32 erg. 

We note that while we can detect cavities everywhere 
in the cloud, the opacity of the 12 CO line makes it dif¬ 
ficult for lower velocity or weak outflows to be detected. 
For example, some HH objects are not seen with high 
velocity 12 CO counterparts. Als o the area of the sky 
covered bv lKnee fc SandeH <1200011 is smaller than of our 
13 CO velocity cube and didn’t include sources such as 
IRAS 1 because of this. It is possible that the total mo¬ 
mentum budget f rom outflows is somewh at higher than 
that estimated bv lKnee fc Sandelll <1200011 . 

We now compare the energy flux from outflows into 
the cloud to that being dissipated in the cloud. Since 
turb ulent energ y decays in roughly a crossing time 
llElmegreen fo Scald2004 and references their-in), the en- 
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ergy dissipated by a turbulent medium can be estimated 
as 


L 


turb 


M c i our ivf urb 


l 


eddy 


(7) 


where l e ddy is the size of the largest turbulent eddies. 
Inserting M c ioud = 520 Mm for the mass of the cloud (us¬ 
ing that measured by iRidge et al.l 120031 but correcting 
for distance), Vturb ~ lkm/s and l e ddy ~ 0.3pc we esti¬ 
mate a total energy dissipation rate of L turb ~ 10 33 erg/s. 
We insert l e ddy ~ 0.3 pc corresponding to the approxi¬ 
mate size of the filaments on the cloud’s southern side. 
The estimated energy dissipation rate is close to that 
available from the outflows seen in 12 CO and that in¬ 
ferred from the cavities. This confirms previous stud¬ 
ies that have found that outflows could be responsi¬ 
ble for drivi ng a significant fraction of molecular cloud 
turbulence IBallv et all 1 199611 IBallv fc R emurthl l2001t 

iKnee fe Sandellll2000HWarin et al.lll996fl . However, here 
we have provided an additional intermediate step be¬ 
tween the driving of the outflow and the excitation of 
the turbulence, the slow expansion of multiple cavities 
which have been opened by now extinct YSO outflows. 

There is an additional test we can carry out. Previ¬ 
ously we estimated the volume filling factor of the cav¬ 
ities ~ 10%. We also have predicted how the cavities 
evolve in time. We ask: what would the cavity expan¬ 
sion velocity be when it reaches a size large enough to 
overlap a nearby cavity? When do random turbulent mo¬ 
tions in the cloud begin to destroy the cavities? It may 
not be a coincidence that the cavity expansion velocities 
~ lkm/s are similar to the rms velocity dispersion of the 
cloud. From Equation^we find that for a given momen¬ 
tum impulse and ambient density, the expansion velocity 
is only weakly dependent on the cavity size, v exp oc R 1 / 3 . 
Since the cavity expansion velocity is only weakly depen¬ 
dent on the cavity size, the expansion velocity would be 
similar for larger cavities. The cavities, when they be¬ 
come overlapping, would have expansion velocities near 
1 km/s which is approximately the same as the cloud rms 
velocity or velocity dispersion. Since this velocity is sim¬ 
ilar to the random component of the cloud, the cavities 
will be destroyed by the turbulence of the cloud, which 
they in fact are contributing to. The high volume filling 
factor of these cavities, and similar expansion velocity to 
the rms velocity dispersion of the cloud also suggests that 
they are a significant contributor to the cloud turbulence. 


5. SUMMARY AND DISCUSSION 

In this paper we began by studying the structure of 
the integrated properties of the 13 CO velocity cube in 
comparison to the distribution of young stellar objects 
seen at near-infrared wavelengths. The velocity disper¬ 
sion of the molecular cloud is between 1 and 1.6 km/s 
rms and varies little across the cloud and does not seem 
to be closely related to the distribution of young stellar 
objects. The mean velocity does show structure, with 
the southern half of the cloud at a mean velocity that is 
about 0.3km/s below that of the northern half. Little ev¬ 
idence for high velocity gas is seen in 13 CO velocity cube 
near the large concentration of young stellar objects or 
outflows identified from other studies. 

More careful inspection of the 13 CO velocity cube re¬ 
veals what appear to be numerous cavities and shells. 


These depressions are likely to be real because they are 
seen in multiple channels, are present in both channel 
maps and position velocity diagrams, and some are also 
seen in the lower signal to noise C 18 0 velocity cube. 
Rims are limb brightened in individual channel maps 
but also tend to contain brighter gas at higher and lower 
velocities than the mean. The features in the position 
velocity diagrams suggest that the cavities may not be 
spherical. They could be expanding ovate or cylindrical 
structures. Examples of cavities likely to be at different 
orientations are seen. Typical cavity sizes are 0.1-0.2 pc 
in diameter and have velocity widths 1-3 km/s. From 
the channel maps we estimate the density contrast be¬ 
tween gas outside and inside the cavity is low, only a 
factor of 2-4. However, the true density contrast could 
be higher if the edges of the cavities are turbulent and so 
contain gas at a range of velocities. We have identified 
22 cavities in a 1 x 1 pc region on the sky. We estimate 
their volume filling factor to be ~ 10%. They permeate 
the molecular cloud. 

We find three cavities with higher velocity 13 CO emis¬ 
sion (above 2 km/s from the cloud mean), suggesting 
that some cavities could be younger than others. Some 
cavities appear to be connected to others. Tubular cavi¬ 
ties are oriented north-south on the southern part of the 
cloud and east-west on the eastern side. These cavities 
could have been driven by previously active outflows lo¬ 
cated near the center of the cloud where there are million 
year old star clusters. 

We have had difficulty matching cavities with young 
stellar objects that could have been responsible for open¬ 
ing the cavity. Cavities are not seen near HH objects or 
objects identified as outflows from previous surveys. Nei¬ 
ther do the cavities appear to be clearly associated with 
million year old pre-main sequence objects which would 
be present in the 2MASS Ks-band images. We estimate 
that the timescale for a cavity to fill in at the sound speed 
or at the random turbulent cloud velocity would not be 
long, less than a million years, thus the cavities must 
have been produced within the past million years. If 
cavities are relics associated with previous outflow activ¬ 
ity, then they could be located distant from their driving 
source. Also their source stars could have drifted from 
their original location. A star moving lkm/s drifts out 
of a O.lpc sized cavity associated with its formation in 
only 10 5 years. 

From scaling laws, we find that expanding bubbles 
driven by a jet impulse can match the rough physical 
properties of the cavities; i.e., their expansion velocities 
and sizes. The cavities are most likely relics associated 
with previous outflow activity, as previously suggested 
bv IBallv et alJ lll999f) in the case of the Circinus cloud. 
We use the dimensional scaling relations to estimate the 
age and jet/wind impulse total momentum required to 
make a cavity. This total momentum is about 1 Mg,km/s 
and the estimated cavity age is a few hundred thousand 
years. The age estimate is approximately consistent with 
the number of cavities, ~ 20, and the number of active 
outflows in the cloud seen at higher velocities in 12 CO. 

From a comparison between the momentum flux esti¬ 
mated from active outflows measured from higher veloc- 
itv (greater than a few km/s from the mean) 12 CO by 
(Knee fc Sandell ; 200 0; and the requirements to open the 
cavities we estimate that most of momentum from out- 
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flows is fed back into the cloud. This total mechanical 
energy is a significant fraction of that required to power 
the cloud turbulence. However, here we have provided an 
additional intermediate step between the driving of the 
outflow and the excitation of the turbulence; the slow ex¬ 
pansion of multiple cavities which have been opened by 
now extinct YSO outflows. The similarity between the 
cavity expansion velocity widths and the cloud velocity 
dispersion also suggests that the cavities could form an 
intermediate step between the outflows and the excita¬ 
tion of the cloud turbulence. 

Our estimates of the momentum requirements of these 
cavities are based on rough order of magnitude scaling 
arguments. Our preliminary simulations have been re¬ 
stricted to jet impulses, rather than jets followed by the 
slow opening of a cavity. Better models of a longer 
interaction between jets and winds from young stellar 
objects with the cloud may allow better constraints on 
the properties of the cavities and on how they impact 
the evolution of the molecular cloud. Future studies 
can explore in more detail the types of wind blow mod¬ 
els which could account for these cavities, and by fit¬ 
ting observations (channel maps) of individual cavities 
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TABLE 1 

Cavities identified in 13 CO(1-0) 


Name 

RA(J2000) 

DEC(J2000) 

Velocity 

Vel. Range 

Cl 

3:29:14 

+31:09:00 

5.5 

5.2-7.0 

C2 

3:28:29 

+31:03:20 

6.0 

-6.5 

C3 

3:29:03 

+31:05:50 

6.1 

5.5-6.5 

C4 

3:29:14 

+31:23:00 

6.3 

-6.8 

C5 

3:29:14 

+31:15:20 

6.3 

4.5-7.0 

C6 

3:28:19 

+31:01:50 

6.8 

6.5-7.5 

C7 

3:29:29 

+31:20:10 

6.9 

6.5-8.0 

C8 

3:28:32 

+31:16:10 

7.0 

-7.5 

C9 

3:29:16 

+31:06:00 

7.3 

6.5-8.0 

CIO 

3:28:58 

+31:18:00 

7.3 

5.5-10.5 

Cll 

3:28:50 

+31:07:00 

7.3 

6.5-7.8 

C12 

3:28:51 

+31:02:00 

7.5 

6.5-7.8 

C13 

3:28:37 

+31:11:30 

7.6 

7.0-8.5 

C14 

3:29:09 

+31:12:00 

7.7 

7.0-8.3 

C15 

3:29:45 

+31:12:10 

8.0 

7.0-9.5 

C16 

3:29:25 

+31:26:30 

8.0 

7.5- 

C17 

3:29:10 

+31:24:50 

8.1 

7.5-8.5 

C18 

3:29:22 

+31:15:30 

8.3 

7.8-9.0 

C19 

3:29:51 

+31:18:10 

9.1 

8.5-11.0 

C20 

3:28:52 

+31:15:50 

9.1 

8.5- 

C21 

3:29:05 

+31:14:20 

9.1 

8.0- 

C22 

3:29:09 

+31:24:10 

9.3 

8.5- 


Note. — Cavities are labeled in 13 CO(1-0) channel maps and position velocity diagrams and channel maps shown in Figures 14151 and 
[HI Central cavity RA and DEC are given with EPOCH J2000. The central velocity is given with respect to the local standard of rest and 
is in km/s. The rightmost column shows the approximate mininum and maximum velocities (in km/s) of the channels in which the cavity 
is seen. 
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Fig. 1.— Mean velocity of the 13 CO(1-0) velocity cube contours fsee lRidge et al.l2003D overlayed on the Ks band 2MASS image shown 
as grayscale. Contour spacing is 0.3km/s and given with respect to the local standard of rest. Red contours are at 8.2 km/s and higher. 
The yellow contour is at 7.9 km/s. Green counters are at 7.3 and 7.6 km/s. The blue contours are at 7.0 km/s and lower. The bright region 
just to the left a nd below the center of the frame contains the star SVS13 which has been used as a reference point by previous works (e.g., 
IKnee fc Sandelill200flr) . Also shown for reference are the B stars BD+30547 and BD+30549. BD+30549 illuminates the reflection nebula. 
The region south of SVS-13 is moving at a higher mean velocity (by about 0.5km/s) than the northern region. The angular resolution of 
the 13 CO data is 47”. Coordinates are given with respect to epoch J2000. 
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Fig. 2.— Velocity dispersion (rms) about the mean from the 13 CO(1-0) velocity cube are overlayed on the integrated 13 CO intensity 
(grayscale). Contour spacing is 0.3 (km/s). The red contours are at 1.5 and 1.8 km/s. The yellow contour is at 1.2 km/s. The green contour 
is at 0.9 km/s. Much of the structure in the velocity dispersion map is due to the low signal to noise at the cloud outer edges. While the 
cloud shows signifi cant variations i n the mean velocity, little variation in the vel ocity disp ersion is seen across the cloud. For comparison, 
IRAS sources from l.len ning^t^rlj ar e show n as red stars. HH objects from lReipurthl (1199911 are shown as orange diamonds. Compact 

submillimeter sources from lKnee fc Sandelll i20Q?1) are shown as blue triangles. The B stars BD+30547 and BD+30549, and infrared source 
SVS13 are shown as dark blue crosses. Two depressions in the molecular gas distribution are seen, one just to the north west of SVS13 and 
the other 7 areminutes north and slightly east of SVS13. These depressions are not coincident with bright infrared sources, HH objects or 
early type stars. 









Right Ascension 


Fig. 3. — Integrated 13 CO intensity (as contours) overlayed on top of the Ks band 2MASS image (as grayscale). Contours are shown 2 K 
km/s apart with the lowest and highest contours at 2 and 14 K km/s. The 2MASS image shows T-Tauri cluster stars as well as background 
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Fig. 4. — Channel maps are shown for a subsample of the velocity cube. The velocity spacing between the channel maps displayed is 
4 channels or 0.53 km/s. Each channel has a velocity width of 0.133 km/s. The velocity of the channel shown in each panel is shown on 
the top right hand side of the image. Multiple shells and cavities are evident in the individual channel maps. Depressions in the intensity 
identified in both position velocity diagrams and individual channel maps are labeled here on the channel maps. The properties of these 
cavities are listed in Table 1. 
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Fig. 5.— Position velocity maps. Each map is located at a different declination in the 13 CO(1-0) velocity cube. This DEC is shown 
on the top right hand side of each panel. Velocity is shown as the vertical axis, whereas the horizontal axes are RA(J2000). The contour 
spacing is 0.5K with the lowest contour at 0.5K. Depressions in the intensity identified in both position velocity diagrams and individual 
channel maps are labeled here on the channel maps. The properties of these cavities are listed in Table 1. 
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Fig. 6.— Position velocity maps. Each map is located at a different right ascension in the 13 CO(1-0) velocity cube. This RA is shown 
on the top right hand side of each panel. Velocity is shown as the vertical axis, whereas the horizontal axes are DEC(J2000). The contour 
spacing is 0.5K with the lowest contour at 0.5K. 
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Fig. 7. — Position velocity maps from the C 18 0(1-0) velocity cube. Similar to Figure 0 The contour spacing is 0.2K with the lowest 
contour at 0.2K. Cavities in the denser regions of the cloud identified from the 13 CO(1-0) velocity cube are also evident in the C ls O 
position velocity maps. 
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Fig. 8. — Positions of cavities compared to the locations of other sources. IRAS sources from IJennings e~ :& q cap are sh o wn as red 
stars. HH objects from iReiourthl J1999T) are shown as orange diamonds. Compact submillimeter sources from iKnee XSandelll d200CT) are 
shown as blue triangles. The cavities are shown as green circles of radius 1.5 arcminute. The grayscale background is the Ks band 2MASS 
image and shows T-Tauri cluster stars as well as background stars. The B stars BD+30547 and BD+30549, and infrared source SVS13 are 
shown as dark blue crosses. 











